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Abstract 

We explore the possibility that the excess neutral and charged current events 
seen by the HI and ZEUS Collaborations at HERA at large x and could 
arise from a hard charm component of the nucleon. While the symmetric 
intrinsic charm hypothesis is unable to account for the HERA anomaly, a 
non-symmetric charm distribution generated non-perturbatively, for which c 
is much harder than c, can produce significant enhancement of cross sections 
at the HERA kinematics. 

PACS numbers: 13.60.Hb, 14.65.Dw, 12.38.Qk 
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The recent observation by the HI and ZEUS |0] Collaborations at HERA of an excess 
of events at large x and in e~^p neutral current (NC) and charged current (CC) deep- 
inelastic scattering (DIS) has prompted numerous speculations about whether one has seen 
evidence for physics beyond the standard model. Most explanations have concentrated on 
leptoquarks or contact interactions as the source for the anomalously large cross sections 
— for a recent overview of the theoretical discussions see Ref. 0. Using a simple, non- 
perturbative model for the intrinsic charm component of the nucleon wave function, we find 
a hitherto unnoticed asymmetry, with the c component significantly harder than c, which 
has important consequences for the HERA NC and CC events. Confirmation of such an 
asymmetry would challenge our understanding of the non-perturbative structure of hadrons 
within QCD. 

The ZEUS experiment collected only NC data, while HI observed both NC and CC DIS 
events. The evidence for the excess was found in data with ^ 10000 GeV^, and with an 
invariant mass M = s/xs of the virtual photon — parton system above M ~ 150 GeV, where 
X is the Bjorken variable, and s ~ (300 GeV)^ is the e^p center of mass energy squared. 
For > 35000 GeV^ two events were observed in the ZEUS data while ~ 0.15 were 
expected from the standard DIS model, and for x > 0.55 and y = M'^/Q'^ > 0.25 four events 
were found with ~ 0.9 anticipated. For > 15000 GeV^ the HI data sample contained 12 
NC (4 CC) events, where ~ 4.7 (1.8) were expected. On the face of it, this would amount 
to a probability of ^ 1% that the excess events arise from statistical fluctuations within the 
standard model. On the other hand, it has been argued 0| that the HI and ZEUS data sets 
are at present as incompatible with each other as each is with the standard DIS model, so 
that more data analysis would be needed before definitive conclusions could be reached. 

Given the importance of the potential discovery of new physics, if the excess persists in 
future HERA data it will be essential that one exhausts all possibilities that the phenomenon 
can be accommodated within standard model parameters. Recently, Kuhlmann, Lai and 
Tung have suggested that an increased valence u quark distribution at large x could 
enhance the deep-inelastic cross sections appreciably. The suggested modification to the 
valence u distribution, uv{x) — >• uv{x) + Su{x), where 

6u{x) =0.02 {l-x)°-\ (1) 

is rather hard, as illustrated in Fig.l. Kuhlmann et al. |^ observed that through the QCD 
evolution feed-down effect, evolving this component from ~ 25 GeV^ to ~ 40000 GeV^ 
increases the standard model cross section significantly beyond x = 0.75, without sacrificing 
the quality of the global fits to existing data at lower energies. 

While a distribution of type ([l|) can lead to better agreement with the HERA data, the 
physical origin of such a hard component of the u quark distribution is difficult to motivate, 
given the small mass of the u quark. In addition, it was argued in a recent reanalysis of 
SLAG data for x > 0.7 that the proposed modification would overestimate the large-x 
data by some two orders of magnitude. 

Another possibility raised in Ref. is that the charm quark distribution might be 
enhanced at large x. Because the couplings of u and c quarks to electroweak bosons are 
identical, the effect of a modified c distribution would be the same as that of a modified 
u. In addition, one could more easily imagine that the charm quark distribution could be 
rather hard, owing to the large c quark mass. Such a component would of course have 
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to be generated non-perturbatively, since the charm quark distribution arising from gluon 
bremsstrahlung alone looks like a typical soft sea distribution, namely peaks at x — and 
is negligible beyond x ~ 0.4, Fig.l. 

The effect of a non-perturbative, or intrinsic, charm component on the DIS cross sections 
at HERA kinematics was recently investigated by Gunion and Vogt |^ within a model of the 
5-quark component of the nucleon wave function on the light-cone [§. Following Brodsky et 
al. [^], the wave function was assumed to be inversely proportional to the light-cone energy 
difference between the nucleon ground state and the 5-quark excited state. The resulting 
a;-dependence of the inclusive c quark distribution in the minimal model of |^ was given by 

5(^f^)c(x) = ^1^^ -x){l + lOx + x^) - 6x(l + x) log l/x) , (2) 

with normalization fixed to 1% [0]. As seen in Fig.l, such a distribution peaks at x ~ 0.2, 
and is negligible beyond x ~ 0.7. The anti-charm distribution is assumed to be equivalent 
to the charm distribution in this model, 5'^^'"^c{x) = 6^^'"^c{x). After evolving to the HERA 
kinematics, this intrinsic charm distribution, while considerably harder than that generated 
through pQCD alone, is still too soft to account for the excess HERA events [0. 

What one needs, therefore, is a substantially harder distribution which has significantly 
more strength above x ~ 0.6 than in (|^). In this paper we investigate whether one can 
indeed obtain a sizable enhancement of the number of events at large x and from a hard 
charm distribution, and whether such a distribution can be motivated from any physical 
model. 

As an alternative to the intrinsic charm picture of Refs. the charmed sea was taken 
in Refs. [§,0] to arise from the quantum fluctuation of the nucleon to a virtual + Ac 
configuration. The nucleon charm radius ^ and the charm quark distribution |10] were both 



estimated in this framework. The meson cloud model for the long-range structure of the 
nucleon has been developed in Refs. |TlHl3i to describe various flavor symmetry breaking 
phenomena observed in DIS and related experiments. It offers a natural explanation of 
the d excess in the proton over u ||14| in terms of a pion cloud, which itself is a necessary 



ingredient of the nucleon by chiral symmetry. It also provides an intuitive framework to 
study the strangeness content of the nucleon, through the presence of the kaon cloud of 
the nucleon [^,0. Whether the same philosophy can be justified for a cloud of heavy 
charmed mesons and baryons around the nucleon is rather more questionable given the 
large mass of the fluctuation. Nevertheless, to a crude approximation, we may take the 
meson cloud framework as an indicator of the possible order of magnitude and shape of 
the non-perturbative charm distribution. Furthermore, a natural prediction of this model is 
non-symmetric c and c distributions. 

In the meson cloud model, the distribution of charm quarks in the nucleon on the light 
cone at some low hadronic scale is written in convolution form: 

S-cix) = [ ^fn/Niz) (^) , 5c(x) = £ ^/a./^(z) (^) , (3) 

where z is the fraction of the nucleon's light-cone momentum carried by the D~ meson or 
Ac. The light cone (or infinite momentum frame) distribution of D~ mesons in the nucleon 
is given by: 
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iojNKz) - ak^— — — , (4j 



I? 



and can be shown to be related to the hght-cone distribution of Ac baryons, fA^/Ni^), by 

fh.lN{z) = fD/N{l - Z). 

In Eq.@ the function g describes the extended nature of the DKcN vertex, with the 
momentum dependence parameterized by ;z) = (A^ + M^)/(A^ + S£)Ac); where the 

DAc center of mass energy squared is given by sdk^ = (^1 + M]^)/ z + {k^ + M|^)/(l — z), 
and go is the DAcN couphng constant at the pole, sda^ = ^n- We expect to be similar 
to the TcNN coupling constant. 

Because of the large mass of the c quark, one can approximate the c distribution in the 



D meson |T0[ and the c distribution in the A+ by: 

c^'{x) ^6{x-l), c^^(x) ^ 6{x-2/3), (5) 



which then gives: 
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5c{x) ^ fo/Nix), Sc{x) ^ -/A^/iv(3x/2). (6) 



The resulting Sc and Sc distributions are shown in Fig.l, calculated for an ultraviolet 
cutoff of A ^ 2.2 GeV, which gives Jq dx6c{x) = dx6c{x) ~ 1%. For a probability of 0.5% 
one would need a smaller cutoff, A f» 1.7 GeV. Quite interestingly, the shape of the c quark 
distributions is quite similar to that in the intrinsic charm model of Refs. 0,^]. However, 
as mentioned above, the model of assumes identical shapes for the non-perturbative c 
and c distributions, while the meson cloud gives a significantly harder c distribution. 

Calculation of the NC and CC cross sections requires parton distributions for all flavors. 



For this we use a recent parameterization of global data from the CTEQ Collaboration [16 



Expressions for the differential NC and CC cross sections d'^a/dxdQ'^ in the standard model 



can be found in Refs. and [0. In Fig. 2 we show the ratios of the modified to standard 
DIS model NC cross sections, with a = d'^a/dxdQ'^, and a + 6a represents the cross section 
calculated with the modified distributions. The result with the modified u distribution 
rises sharply above M ~ 200 GeV. The effect is rather similar if one includes the non- 
perturbative 6c and 6c distributions from the meson cloud model. On the other hand, with 
the somewhat softer, intrinsic charm distribution of the enhancement is rather modest, 
and ^ 10% over the whole range of M. 

Also shown in Fig. 2 is the effect of additional contributions to the light flavor distributions 
{u, u, d and d) which one would obtain from a pion cloud of the nucleon [|ll|. Although these 



are considerably larger in magnitude than the 6c or 6c distributions, because they appear 
at small x (~ 0.1) their effect is to yield only a very small enhancement of the cross section 
ratio. From this figure one can conclude that the only realistic candidates for a significant 
enhancement of the cross section at large M are the modified valence u distribution from 
Ref. 0, and the hard charm distributions in Eq.(^. 

To compare with the HI and ZEUS data one needs to evolve the cross sections to the 
typical scales relevant for HERA. Since the excess of events is seen at quite large Q^, 10000 ^ 
< 40000 GeV^, we evolve the 6u and 6c & 6c distributions to an average value of 
= 20000 GeV^. The result is shown in Fig. 3, where for M > 250 GeV especially one 
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has an enormous enhancement of the NC cross section. The enhancement is quahtatively 
similar for both the 0.5% and 1% 5c components. 

For scattering via the CC, the effect of the additional contributions to the parton distri- 
butions is shown in Fig.4 for = 20000 GeV^. Since the boson is not sensitive to the u 
quark in the proton, the 6u modification has no effect on the e~^p cross section. In contrast, 



as noted by Babu et al. ||T8[, the effect of an additional non-perturbative 5c contribution 
is an even larger enhancement of the CC cross section than the NC cross section. With a 
0.5% (1%) intrinsic charm component the CC cross section increases by a factor ~ 2 (3) for 
200 ^ M ^ 250 GeV, which is similar to the excess observed by HI []T| in this region. 

Having seen that one can indeed modify the existing parton distributions at large x to 
produce large enhancements of the cross sections at large M and Q^, the key question to 
address next is whether such modified distributions are compatible with existing data at 
lower Q^. A recent analysis of SLAC data at 0.7 < a; < 0.97 in and near the resonance 
region suggests that the proposed modification of the u quark distribution in Ref. 0] would 
overestimate the large-x data by two orders of magnitude. (Though one should keep in 
mind that the authors of Ref. assume the validity of Bloom-Gilman duality even at very 
large x values close to 1.) On the other hand, the values of corresponding to the x 
and bins are too low for the SLAC data to be sensitive to the charm component of the 
nucleon wave function. Data are, however, available from the BCDMS Collaboration in the 



region x > 0.6 which are above charm threshold ||T9[. The effect of adding the intrinsic 



6c distribution calculated from Eq.(|) to the data on the deuteron structure function (as 
parameterized in Ref. pOf) is illustrated in Fig.5. With the addition of the 0.5% contribution 
there does not seem to be any inconsistency with the data [0 and, indeed, the agreement is 
slightly improved at x = 0.75. On the other hand, the 1% case may be a little too high for 
comfort. Earlier data from the EMC PT| on the charm structure function extracted from 



charmed particle production in /i-nucleus scattering also exist, but cover only the region 
below X ~ 0.3. 

One can conclude, therefore, that a non-perturbatively generated hard charm distribution 
can significantly enhance the DIS cross sections at the large x and values attainable at 
HERA, and still be compatible with data at lower Q^. We should make clear, however, that 
the Sc and 6c quark distributions calculated in the meson cloud model should not be viewed 
as quantitative predictions, but rather as a physical motivation for a hard non-perturbative 
charm component of the nucleon. In particular, one can vary model parameters to modify 
the detailed shape of the distributions. For example, using a smaller power of 1/(A^ + s_daJ 
in the hadronic vertex function gihPi^z) would result in a 6c distribution with even more 
strength at a; > 0.75 than that in Fig.l, which would further enhance the cross sections 
for M ~ 200 — 250 GeV. Nevertheless, the fact that the model predicts a large asymmetry 
between c and c distributions is a robust, and to our knowledge unique, feature of the meson 
cloud model. 

Another possibility which could lead to additional enhancement of the cross sections at 
large x would be a non-perturbative h quark component , although the magnitude of this 
would presumably be rather small compared with the c. In addition, the recent reanalysis 
of the SLAC deuteron data in Ref. |22| suggested that the valence d/u ratio does not tend 
to zero as x — 1, as assumed in global fits to data |T6|, but rather is consistent with the 
expectation of perturbative QCD, namely d/u — » 1/5 Figure 6 shows the ratio of the 
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modified d quark distribution, d + Sd, extracted from tlie deuteron and proton data in Ref. 



H at ~ 10 GeV^ and evolved to = 10000 GeV^, to tliat used in the standard fits. 



The correction Sd is taken from the parameterization in Ref. . The effect of the modified 
d quark is not large in NC events, primarily due to the factor 4 charge suppression compared 
with u, but can be sizable for CC events, comparable to that of the 0.5% 6c scenario — as 
shown by the dashed curve in Fig.4. The fact that all standard sets of parton distributions 
assume that d/u^Oasx^l means that there is a possible source of systematic error in 
the modeling of "background" events which should be accounted for. 

The task of disentangling the origin of the HERA anomaly in future can be approached 
from several directions. Tagging charm final states in J/ip production could provide infor- 
mation on whether the large x enhancement is due to charm or other flavors. A thorough 
comparison of NC and CC cross sections, for both e~^p and e~p collisions, would also enable 
one to determine the flavor combinations responsible for the cross section enhancement 



One can also construct new distributions, such as those suggested in Ref. |^5|], which are 
insensitive to parton distributions, but sensitive to modifications of parton level cross sec- 
tions, and which could therefore distinguish between explanations of the anomaly in terms 
of modifications of parton distributions and those invoking new physics. 

Needless to say, we eagerly await further results from the HI and ZEUS Collaborations. 
Quite apart from the excitement over the possible discovery of physics beyond the standard 
model, the issue of whether or not the nucleon contains a significant component of intrinsic 
charm is extremely important. In particular, any experimental evidence supporting the 
suggestion that the intrinsic charm could have a strong asymmetry would mean a revision of 
current wisdom, and would undoubtedly lead us to a deeper understanding of the structure 
of hadronic systems. Whether HI and ZEUS provide evidence for new physics beyond the 
standard model or map out hitherto unknown features of the non-perturbative structure of 
the nucleon, the results will be of great interest. 
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FIG. 1. Modification 6u of the valence u distribution from Ref. compared with the hght-cone 
intrinsic charm distribution, of Refs. [^Q, and the non-perturbative 5c and 6c distributions 

in the meson cloud model (the latter normalized to 1%). Also shown is the purely perturbative 
contribution at =4 GeV^. 
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FIG. 2. Ratio of modified to standard DIS model NC cross sections at the charm threshold, with 
the modifications arising from the additional u quark component (dashed), 1% non-perturbative 
6c and 6c distributions from the meson cloud model (solid), and the intrinsic charm model of Refs. 
I^JH (dot-dashed). Also shown is the effect of the meson cloud contributions to the light sea quarks 
111 (dotted). 
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NC Q^ = 20000GeV- 




modifications are from tlie extra 5u component [||, and from additional hard non-perturbative 6c 
and 5c components, with 0.5% and 1% normahzation. 
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FIG. 4. Ratio of modified to standard DIS model CC cross sections at = 20000 GeV^, with 
the modifications arising from 0.5% and 1% additional 6c distributions, as well as a modified d 
quark distribution at large x (see Fig.6). 
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FIG. 5. Structure function of the deuteron, F2D, with 0.5% (sohd) and 1% (dashed) additions 
of non-perturbative charm to the global fit (dotted) from Ref. . 
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FIG. 6. Ratio of the modified d quark distribution, d + 6d, extracted from deuteron data in 
Ref. H at Q2 ~ 10 GeV^ and evolved to = 10000 GeV^, to that used in standard, global fits 
|16|, which assume d/u ^ as x — > 1. 
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